SEMICLASSICAL LIMIT FOR THE SCHRODINGER EQUATION 
WITH A SHORT SCALE PERIODIC POTENTIAL 



FRANK HOVERMANN, HERBERT SPOHN, STEFAN TEUFEL 



Abstract. We consider the dynamics generated by the Schrodinger operator 
H = - |A + V{x) + W(ex), where V is a lattice periodic potential and W an 
external potential which varies slowly on the scale set by the lattice spacing. 
We prove that in the limit e — > the time dependent position operator and, 
more generally, semiclassical observables converge strongly to a limit which is 
determined by the semiclassical dynamics. 



1. Introduction 

A basic problem of solid state physics is to understand the motion of electrons 
in the periodic potential which is generated by the ionic cores. While this problem 
is quantum mechanical, many electronic properties of solids can be understood 
already in the semiclassical approximation H, 16, One argues that if the wave 
packet spreads over many lattice spacings, the kinetic energy (/ifc)^/2m is modified 
to the n-th band energy £'„(fc). Otherwise the electron responds to external fields, 
Eex, Bex, as in the case of vanishing periodic potential. Thus the semiclassical 
equations of motion are 

f = Vnik) = VkEn{k) 

(1.1) 

hk = e(Sex(r) + Vn{k) A B^^{r)) , 

where r is the position and k the quasimomentum of the electron. Note that there 
is a semiclassical evolution for each band separately. 

The goal of our paper is to understand on a mathematical level how these semi- 
classical equations arise from the underlying Schrodinger equation. We consider 
only the case where i?cx = 0. 

The setup is rather obvious. We start from the Schrodinger equation 

(1.2) '^'^^^^ 
with Hamiltonian 

(1.3) H = -]^/:\ + V{x) +W{ex). 



The electron moves in R'' and the solution to (1.2) defines the unitary time evolution 
U^(t)'ip{x) = e~^*^il}{x) = ip{x,t) in L'^{R'^). We have chosen units such that h = 1 
and the mass of the particle m = 1. ^(a;) is a periodic potential with average 
lattice spacing a. The precise conditions on V will be spelled out in the following 
section, where we also describe the direct fiber integral decomposition for periodic 
Schrodinger operators. The lattice spacing a defines the microscopic spatial scale. 
W{ex) is an external electrostatic potential with dimensionless scale parameter £, 



2 



FRANK HOVERMANN, HERBERT SPOHN, STEFAN TEUFEL 



£ ^ 1, which means that W is slowly varying on the scale of the lattice. For real 
metals the condition of slow variation is satisfied even for the strongest external 
electrostatic fields available, cf. Chapter 13. 

The external forces due to W are of order e and therefore have to act over 
a time of order to produce finite changes, which defines the macroscopic time 



scale. We will mostly work in the microscopic coordinates {x,t) of (1.2). For sake of 
comparison we note that the macroscopic space-time scale (x' , t') is defined through 
X = e~^x' and t = e^^t' . With this scale change Eqs. (1.2), (|1.3|) read 



d 



(1.4) 



H 



-e\A' + V{x'/e) + W{x') 



with initial conditions ip^{x') = £~'^^'^'\\){x' je). If ^ = 0, Eq. (F4) is the usual 
semiclassical limit with e set equal to ti. Thus our problem is to understand how an 
additional periodic, but rapi dly os cilla ting potential modifies the standard picture. 

The two scale problem (L2), (O) can be attacked along several routes. A 
first choice would be time dependent WKB H, |[ |l^, |l^. In the limit e ^ 0, 
for each energy band separately, one obtains a Hamilton- Jacobi equation for the 
phase and a transport equation for the amplitude of the wave function ?/'(x,t). As 
a main draw-back of this method, generically, the solution to the Hamilton- Jacobi 
equation develops singularities after some finite macroscopic time. If = 0, it is 
well understood how to go beyond such caustics by introducing new coordinates 
on the Lagrangian manifold. For ( |l.2|) , (1.3) a corresponding program has not yet 
been attempted. The results |5[ |6|, |10|, |12|] are valid only over a finite macroscopic 
time span with a duration depending on the initial wave function. 

Another variant is to establish the semiclassical limit through the convergence 
of Wigner functions. In our context one defines a band Wigner function W^{r, k,t) 
depending on the band index n and as a function of the position and quasimo- 
mentum. One then wants to prove that in the limit e — )■ W^{t) converges to 
Wn{t), which is the initial band Wigner function VF„(0) evolved according to the 
semiclassical flow (1.1). Such a result is established in ^, |8| for the case of zero 
external potential, the general case being left open as a challenging problem. 

A third approach to the semiclassical limit for V — is the strong convergence of 
Heisenberg operators ^, |l^, |2^. We briefly recall its main features. We define, 
as unbounded operators onL?(Mr), 



x{t) :-- 



pit) := e'*"pe-'*", p = -iV^, 

where H is the Hamiltonian in (^^) with V — 0. The goal is to establish the strong 
limit of 

x''{t)ip = ex{e~'^t)ip, p^it)ijj = p{s~'^t)ip 

as e ^ with ip in a suitable domain. In the trivial case of free motion, W = 0, 
this amounts to the strong convergence of x'^{t)ip = {ex +pt)'ip^ p'^{t)tp = ptp^ which 
yields lime^ox^(i) = pt, lime^op'^(i) — P- The general case requires more work 
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22 1 . One obtains the strong limits 



(1.5) 



linia;^(t) = r{p,t) , 

e—fO 



e— >0 



Here r{p,t), u{p,t) are solutions of 

(1.6) f = M, u^-SIW(r) 

with initial conditions — Q, u^) — p. The initial condition — Q reflects that 
IV'P looks like 5{r) on the macroscopic scale, provided that HV'lb = 1- For general 
initial conditions, tq ^ 0, we would have to shift the initial ■0 by e^^rg. 

The strong operator convergence may look slightly abstract, but all the desired 
physical information can be deduced. E.g., the initial '0 defines the momentum 
distribution |?/'(A:)p independent of e and the 5{r) spatial distribution in the limit 



e — > 0. Then, according to (1.5), for small e the position distribution at time t is 
given by 

~ {^J{r{p,tm ^ jmk)\y{rik,t))dk, 

which means that the phase space distribution 6{r)\'i/j{k)\'^ dr dk is transported ac- 
cording to the semiclassical flow ( |r^) . The spatial marginal of this distribution 
at time t is the desired approximation to the true position distribution |'(/;^(x, t)|^. 
\'ip'^{x,t)\'^ may oscillate rapidly on small scales and some averaging, as embodied 
by the test function /, is needed. 

In this paper we investigate the semiclassical limit (^]^) , (^^) through the strong 
convergence of the position operator x'^{t). We will show that, in the limit £ ^ 0, 
x'^{t) is diagonal with respect to the band index and in each band the structure is 



analogous to ( p_.5| ) with p replaced by the quasimomentum k and (1.6) replaced by 
( pTl] ) . More generally we will consider the semiclassical limit of the Weyl quantized 
operators {ex,p), whose classical symbol is periodic in p. 

To give a short outline: In the following section we collect some properties of 
periodic Schrodinger operators. In Section ^ we state our main results, which are 
proved in Sections |[ |[ and ^ respectively. In Section ^ we discuss some impli- 
cations for the position and quasimomentum distributions, and, more generally, 
for the band Wigner functions. The difficulties arising from band crossings are 
explained in Section ^. 



2. Periodic Schrodinger operators 

For the periodic potential V we will need only some rather minimal assumptions, 
which we state as 

Condition (Cpcr)- Let F ~ Z'' be the lattice generated by the basis {71, . . . ,7^}, 
7,; e R*^. Then V{x + 7) = V{x) for all x G M'^, 7 £ F. Furthermore, we assume 
V to be infinitesimally operator bounded with respect to Hq. The last condition is 
satisfied, e.g., ifVE LP{AI), where M is the fundamental domain of T, and p — 2 
for d < 3 and p > d/2 for d > S, respectively. 

(Cpor) will be assumed throughout. 
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We recall the Bloch-Floquet theory for the spectral representation of 
(2.1) Hp,, = ^ _ 

The reciprocal lattice F* is defined as the lattice generated by the dual basis 
{jI,... ,7^} determined by 7^ • 7* — 2TrSij, i,j — 1,... ,d. The fundamental 
domain of T is denoted by M, the one of F* by M*. M* is usually referred to as 
first Brillouin zone. If we identify opposite edges of M, resp. M* , then it becomes 
a flat d-torus denoted by T = W^/V, resp. T* M''/F*. 

Let us introduce the Bloch-Floquet transformation, which should be viewed as 
a discrete Fourier transform, through 

{Ui^){k,x) ^e-^("+'')-'=V(a; + 7), {k,x) e R^rf ^ 

for il) e S{W^). Clearly, 

{Uij){k,x' +^)^[U^){k,x'), 

(2.2) 

{U^){k' + -i*,x) = e"'"-^*(WV)(fc',2:) . 
Therefore it suffices to specify U4> on the set M* x M and, if needed, extend it to 
all of M^d by (U). The linear map U : L2(E^) d 5(M^) ^ H := /®. L^{M)dk, 
with dfc the normalized Lebesgue measure on M*, has norm one and can thus be 
extended to all of L^(M'*) by continuity. U is surjective as can be seen from the 
inverse mapping 

[U-^(j)){x) / e"-''(j){k,x)dk, 
Jm' 

which has norm one. Thus U : L'^{W^) ^ 7i is unitary. 

To transform -ffper under 14, we first note that p = lApU^^ = + fc, with 
Dx = — «Vx- Therefore 

Hp,, := UHp„U-' = / Hp„{k) dk , 
Jm* 

and 

Hp„{k)^^(Dx + kf + V{x), keR". 

Hpe,{k) acts on L'^{M) with fc-independent domain D := H'^{T). ip G D is periodic 
in X. Hp,,{k) is a semi-bounded self-adjoint operator, since by condition (Cper) V is 
infinitesimally operator bounded with respect to — A ||^ . In particular, Hpe,{k) is an 
entire analytic family of type (B) in the sense of Kato for k G C^. Since the resolvent 
of Ho{k) = ^{Dx + k)^ is compact, the resolvent R\{Hpc,{k)) := {Hpe,{k) — A)~^, 
A ^ a{Hpe,{k)), is also compact, and Hpc,{k) has a complete set of (normalized) 
eigenfunctions (pn{k) G iJ^(T), n G N, called Block functions. The corresponding 
eigenvalues En{k), n G N, accumulate at infinity and we enumerate them according 
to their magnitude and multipHcity, Ei{k) < E2(k) < . . . . En{k) is called the n-th 
band function. We note that Hp„{k) = e~" '''*i?per(A; -I- 7*)e" '''' . Therefore i?„(fc) 
is periodic with respect to F*. If i?„_i(fc) < Enik) < En+i{k) for aU k e M* 
(in particular En{k) is nondegenerate), then the n-th band is isolated. In this 
case En and the corresponding projection operator are real analytic functions as a 
consequence of analytic perturbation theory . We denote by X C N the set of 
indices of isolated bands. 
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It will be convenient to have also a notation for the spectral subspaces. Let 
Pn{k) : L'^{M) L?{M) denote the orthogonal projection onto the n-th eigenspace 
of iJper(fc)- Similarly, we set Qn{k) = 1 — P„(A:). Their direct fiber integral is 
denoted by 

Pn= / Pn{k)dk. 

Pn projects onto the n-th band subspace in Ti and P„ = U~^PnU projects onto the 
n-th band subspace in L'^{W^). We have 

{Pn1p){k, ■) = P„(fc)V'(fc, •) = [i^nik), ■lp{k))L2(M)V'n{k, •) 

(2.3) = Mk)Mk,-). 

The coefficient functions '0n G L'^(M*) and are called the Block coefficients in the n- 
th band subspace. For the index set X C N of isolated bands we set Pi = X^nei 

Remark. To have a concise notation, we will use a tilde for operators acting on 
TL. Thus if A is an operator on L'^{W^), then A — UAU^^ . If A has a direct fiber 
decomposition, then A = J^^, A{k)dk with A{k) acting on the fiber L'^{M) ofJi. 

3. Main results 

For the potentials we assume (Cpcr) for V and in addition 

Condition (Cox)- The external potential W G 5(M''). 

To state the semiclassical limit, we first have to explain the classical dynamics 
which will serve as a comparison. For each n € I the classical phase space is 
M'' X T*, where T* ^ R'^/T*. As n-th band Haniiltonian we have 

hn{r, k) = En{k) + W{r), (r, k)eR''xT*, 

and the classical dynamics in the n-th band is governed by 

(3.1) r„ = VkEnikn), k = -VrW{rn). 

Since we want to prove the strong convergence of the position operator, as in the 
case y = 0, we have to lift (^) to operators on H. For this purpose we solve 
( ^ ) with initial condition r„(0) = 0, fc„(0) = k. We denote the solution by 
(r„(i; k), kn(t; k)), regarded as functions of /c G T*. For £ H, we define 

{R{t)ij){k, x) = Y^ rn{t; k)PnMk, x) , 

and analogously, for later use, 

{K{t)i^){k, x)^Yl ^"(^; k)Pni^{k, x) . 

Theorem 3.1. Let the conditions (Cpcr), (Cox) be satisfied. Let 

x^t) ^ eU%-t/e)xU%t/e) . 
Then for every ijj G RanPx D D{\x\) n H'^ , with H'^ the second Sobolev space, 

lim x^ {t)ijj = U-^ R{t)Uip 

e— >0 

strongly. 
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Theorem B.l will be proved in several steps. First we show that in the semiclas- 
sical limit transitions from and to isolated band subspaces are suppressed on the 
level of the unitary groups. We define -ff^iag = PnHPn + QnHQn and U^l^gi'^) ■= 
exp(— iti/^'j^jg). In Section H we will prove 

Theorem 3.2. For any n E T we have 

in B{H^,L^), where is the first Sobolev space. 

The position operator is not diagonal with respect to the n-th band subspace 
and we define its diagonal part by x^jag = PnxPn + QnxQn with the time evolution 

^d'i^gW ^U^~t/e)x^,^^^/s) . 

Our second step is to prove that the off-diagonal part of x"^ (t) vanishes in the limit 
£ 0. 



Theorem 3.3. For n eX 

(3.2) lim{x%t)~x'^.r,^it)) =0 
in B{H^,L^). 

By construction we have [x^dHgi^)^ ^"1 = ^ ^^'^ suffices to study the dynamics 
in the n-th band subspace. This subspace is isomorphic to L^(T*) and, up to errors 
of higher order, x'^^^^{t) can be replaced by xliJ^{t) whose time evolution is governed 
by a Hamiltonian of the form 

H^^^ ^En{k) + W{-ieVk). 

At this stage we can apply the standard machinery of scmiclassics, except that 
formally the roles of position and momentum have been interchanged and the new 
position space is the flat torus rather than M''. 

So far we focused on the position operator, since the electronic density is the 
most accessible quantity experimentally and it corresponds in essence to a suitable 
function of the position. On more general grounds one would like to characterize 
a wider class of semiclassical observables. One further obvious candidate is the 
momentum p. In the Bloch-Floquet basis we have p = k + D^- k is semiclassical, 
as being canonically conjugate to — iVfe: 

Theorem 3.4. Let 

fc^(t) = U^-t/e)U-'^ kUU^t/e). 
Then for every ip € RanPj 

(3.3) lim F(i)V^ = U~^K(t)Ui; 

E— >0 

strongly. 

On the other hand, is unbalanced because there is no extra factor of e. Thus 
p{t/s) has a limit only when averaged over time (compare with Section |q). 



It is relatively easy to see (cf. Section |^) that Theorems 3.1 and imply the 
semiclassical limit also for bounded functions oix^{t) resp. of fc'^(t) (cf. Lemma 3.1). 
Next note that for r*-periodic functions g, g{- +7*) = g{-) for all 7* S F*, we have 
Ug{p)U^^ = g{k) and hence, by the functional calculus for self-adjoint operators. 
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g{p''{i)) = g{k^{t)). Therefore we introduce the set 0(0) C C{W^ x M'^,M) of 
bounded and continuous semiclassical symbols which vanish if the first argument 
approaches infinity and are r*-periodic in their second argument. For a £ 0(0) we 
introduce its Weyl quantization 

(3.4) (aWV)(x) = (^/« (^'^) ^'^''-'^H{y)didy 

as a bounded operator on i^(M'^). The operator corresponding to the symbol 
a{£x,^) will be denoted by a^^^ and we set, as before, 

(3.5) a^^'{t) = U'{~t/e)a:^''W{t/£) . 

Theorem 3.5. Let the conditions (Cpci), (Ccx) be satisfied and a G 0(0). Then 
for every ijj G PjL^ we have 

lim a^''{t)i: U-^a{R(t),K{t))U^p . 
4. Semiclassical distributions 



Theorems 3.1 and 3.5 tell us how the quantum distributions behave in the semi- 
classical limit. Let us first consider the initial V' G PxTi. Its scaled position distri- 
bution is e~'^\ip{x/e)\'^ which converges to S{x) as a measure. The quasimomentum 
distribution X^nei lV'n(fc)P is independent of e. Thus it is natural to chose 

(4.1) p{drdk) = "^S{r)\iPn{k)\^ drdk ^^pn{drdk) 

as initial distribution for the semiclassical flow ( |3.l[ ). We could consider more 
general initial measures at the expense of making ^ itself e-dependent. For example 
the shifted initial measure X]nei^(^ ~ '"o)|V'n(fc)| dr dk is approximated by '4){x — 
£~^ro). Under (3.1) p{dr dk) evolves to p{drdk,t) — X^nei /""("^^ '^^i 0- Each p„ 
satisfies weakly the transport equation 

d 

(4.2) — p„ = ~VEn{k) ■ VrPn + VV{r) ■ VkPn 

at 

with initial condition pn{dr dk,Q) = p„{drdk). We define the position and quasi- 
momentum marginals through 

(4.3) p{dr,t)= / p{drdk,t), p{dk,t) = / p{drdk,t). 

Jm* 

To connect with the quantum evolution we consider the quantum mechanical posi- 
tion distribution 

(4.4) p%dx, t) ^ e-'^\i;{x/e, t/e)\^ dx 



as a probability measure on R'^. From Theorem 3.1 and Lemma B.l we conclude 
that 

(4.5) lim f{x'it))i^ = U-^f{R{t))Ui, 
for / e Coo(K''). In particular, 

(4.6) lim / p%dx, t)f{x) = YiuM, f{x'{m) = {m. f{R{t))m 

E— >0 / e^O 
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and we only have to compute the expression on the right hand side. Using that 

{Uilj){x,k) = ijj,i{k)pn{x,k) 

we have 

(4.7) {U^J{R{t))U^) = J2 f \Mk)\'f{rn{t;k))dk^J2 [ Pn{dr,t)f{r) . 

Thus the positional distribution p'^{dx,t) converges weakly as a measure to the 
incoherent sum ^jjgj /9„((ir, t). 

By the same reasoning, if g is a r*-periodic function, then by Theorem 3.4 and 



Lemma 8.1 



(4.8) lim g{p{t/e))^/j ^ U-^g{K{t))U^ . 

Therefore, if p^{k,t)dk denotes the spectral measure for the quasimomentum op- 
erator at time t/e, we have 



(4.9) limp^(fc,t)dfc = V p„(dfc,t) 

nel 

weakly as measures. 

More generally for ip G L'^ we define the scaled Wigner function by 

(4.10) W'ix,k,t) = ^£-V(e"ix- i7,e-4)^^*(e-ix+ i7,e~4)e*'=-T 

with X eM.'^, k G M*. We think of W as a signed, bounded measure over M'' x M*. 
The Wigner function yields expectations of Weyl quantized operators through 

(4.11) (^,e*^'*/"a^-'e-''-f^'*/^v) = / W%x, k,t)a{x, k) dx dk 



with a r*-periodic in its second argument. From Theorem 3.5 we therefore deduce 
that 

(4.12) VnnW^r,k,t)drdk =p{drdk,t) 



weakly as measures. The limits ( fl.6D and (^^) are the particular cases, where either 
a{x, k) — f{x) or a(x, k) = g(k). 



5. Convergence of the unitary groups 

By definition, the time evolution generated by -ffpcr leaves invariant the band 
subspaces Ran(P„) for all n e N. However, W^{x) = W{ex) does not respect the 
Bloch decomposition and it will induce transitions between different bands. Since 
is of slow variation, we expect such transitions to have a small amplitude as 



stated in Theorem 3.2 
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transforms under U as 
{UW''i;){k,x) = ^e-*(^+'^'-'=iy(£(a; + 7))V(a; + 7) 

= (27r)-'*/2 /■ w{p){Ui}j){k-ep,x)dp 
(5.1) = (271)-'^^'^ J W^p){l(^/j){k^p,x)dp^:{W'Ut{j){k,x), 

where W'^{p) = e~'^W{p/e) and we adopt the quasiperiodic extension (|2.2|). Since 



W e Si'R'^), the integral is well-defined and W = UWU-^ acts on H as 

convolution with in the fiber parameter k. approximates a Dirac delta in 
the limit e — > and the shift in (5.1) becomes the identity operator. 

In the Bloch-Floquet representation the full Hamiltonian (1.3) becomes 

(i/V)(fc, •) - i?pcr(A:)V(fc, •) + {W'4^){k, ■) . 

We expect the diagonal part of to be dominant with the off-diagonal piece 
as a small correction. For such a decomposition it turns out to be convenient to 
fix the index n of an isolated band and to project along P„ and its complement 
Qn 1 — P„. For n GX we define the diagonal part i?d\ag ^ 

^diag = PnHPn + QnHQn , 

and the off-diagonal part of the external potential as 
Then 

H = HS,,^ + = (i7p„. + M/^g) + Cd" ■ 
We note that W^-" and are bounded operators and set 



diag 



To prove Theorem |3 . 2| we start by writing the difference of the two unitary groups 
in the Bloch representation as 

(5.2) il^it/e) - UX,{t/e) ^^le f^'^ W{e-H - s) (e-^C^') U^) ds 
and we have to investigate the operator W^^". By definition, for -ijj E Ti., we have 

{QnW'Pnij){k) = {27r)-^/^ [ W'{p)Qn{k)Pn{k-pMk-p)dp, 



which vanishes strongly i n th e limit e — > 0, since localizes around p ~ Q. To 
control the long times in (5_^) we need uniform convergence of order o(e), however. 
To have a more detailed information on W^^^ we Taylor expand of Pn{k—p) around 
Pn{k), leading, as we will show, to 

(5.3) {QnW'Pn'ikm = -£(2^)-'^/2 / (p)g„ (fc) VfeP„ (fc) ^'(fc -p)dp+ o{e) . 
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Here F'^{p) := W'^{p)^ is the Fourier transform of F'^{x) — {DxW){ex) and we will 
associate to the operator F'^ as in the case of , 

{F'il)){k) = (27r)-''/2 I -p)dp. 



To justify (5.3) we first need to calculate "S/ kPn{k)- 
Lemma 5.1. Let n E T. Then 

VfcP„(fc) = - Qr,ik)RE„ik)iHperikmD, + k)P„ik) 

(5.4) -P„(A:)(D, + fc)i?B„(fc)(i?per(fc))Q„(fc), 

where Rx{H) = (H-X)-^ is the resolvent of H . Thus Pn{-) e C°°{M*; B{L^{M))). 

Proof. Using contour integrals we write 

VfcP„(fc) = I \7kRxiHp,,{k))d\, 

where Cn{k) is a closed rectifiable curve in the complex spectral plane which encircles 
En{k) only. From 

= Vfel = Vfe(i7pe,.(A:)-A)i?A(//pcr(fc)) 

= {D^ + k)Rx{H^,,{k)) + (i/p„.(fc) - A)Vfei?A(i?por(fc)) , 

we infer 

VkRxiHpcAk)) = -RxiHp,Ak))iD.. + k)Rx{Hp,,{k)) . 

Hence we get 

Qn{k)VkPn{k) = Qn{k)VkPn{k){Pn{k)+Qn{k)) 

^ ^ Qn{k)Rx{HpeAk)){D^ + k)Rx{Hp.Ak))Pn{k)d\ 



27^* Jc„(fc) 

= 7T- { ^A(i/por(fc))g„(fc)— dX {D, + k)Pjk) 

^T^l Jc^(k) En[k)~\ 

(5.5) = -RE,^ik){Hpor{k))Qn{k){D, + k)P^{k), 

where the term Qnik)S7 kPn{k)Qn{k) vanishes, since in this case the integrand is 
an analytic function on the whole interior of c„(fc). Note that Pn{k) projects 
onto a subspace of finite energy, on which + fc is bounded. The statement 
about continuity for this term then follows from the continuity of Pn{k), En(k) 
and the assumption that En{k) is isolated from the remainder of the spectrum. An 
analogous computation for Pn{k)V kPn{k) leads to the second term in (5.4). 

Finally, P„(-) £ C°°{M*- B{L^{M))) follows by induction. □ 

From Qn{k) + Pn(k) = 1 we conclude that Qn{k) is differentiable as well and that 
VfcQ„(fc) = -VfcP„(fc). 

Lemma 5.2. Let n El. Then 

W'S - [QnVkPn + PnVkQn) ' F' + o{e) 

m B{n,n), where VkPn := /^^ VfeP„(/c)dfc. 



SEMICLASSICAL LIMIT WITH A SHORT SCALE PERIODIC POTENTIAL 



11 



Proof. We will treat only the QnW^Pn part of W^'^ explicitly, since the argument 
for the second part is analogous. 

Let ip G Tl. By Lemma |3.l| we are allowed to write the following well-defined 
identity, setting Cp = p/\p\j 



W'{p)Pn{k-p)i^{k-p)dp 

W^p)\p\ {\p\~\Pn{k - p) - P„(fc)) + ep ■ VkPnik)) ^P{k - p) dp 
+ P^{k) f W^p)i;{k ~ p) dp ^ [w'{p)p-VkPn{k)^{k-p)dp 



(5.6) = e jw%p)^^l^?^^^^-^-^^^ + ep-V,P^{k)]i^{k-p)dp 
(5.7) 



(27r)'^/2p„(fc)(VK^^)(fc) - e{2iTf'^VkPn{k) ■ (P^^)(fc) . 



If in ( |5.6| ), (5.7) we apply Q„(fc) from the left, the first term of ( |5.7| ) vanishes and 
it remains to show that ( |5.6D , divided by e, tends to zero uniformly for all -0 G ?i. 

We split the integral into two parts. Let i? > be arbitrary, B]^ = {p \ \p\ < R}. 
We start with 



W%p) 



\p\ Pn{k-p)-Pn{k) 



\P\ 



+ ep-VkPn{k) \ i'{k-p)dp 



dk 



< sup sup \\\p\~^Pn{k-p)-Pn{k))+ ep-VkPn{k)\ 

keM' p£Ba 



W'{p) 



M' 



llV^fc -p)IIl2(m) dkdp 



< 



X sup \\\p\-'{Pn{k-p) ~ Pn{k)) + ep-VkPn{k)\\ . 

keM* ,peBR 



Since HF'^H^i does not depend on e and since the difference quotient approaches 
the derivative uniformly on the compact domain M* , the 7i-norm of the first part 
tends to zero uniformly. For the remaining part we have 



rye(M Pnik-p)-Pn{k) 

\p\>R ^ 



M' 



+ Cp ■ VkPnik) ]ip{k -p)dp 



dk 



L2 (Af) 



X sup \\\p\'^{Pn{k-p)-Pn{k))+ep-VkPn{k)\\ , 

keM',peM.'^ 
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which tends to zero uniformly as £ — > 0, since ||i^^||Li(B'=) ~^ for any fixed 
R>0. " □ 



As a consequence of Lemma 5.2 the difference of the two unitary groups in Eq. 
( ^ ) can be written as 

(5.8) = is /o^" U%e-H - s) (Q„VfcP„ + P„VfeQn) • U'^C^is) ds + o(l) . 

We have to estimate the integral without losing one order of e from the integration 
over time. As in the proof in |^ of the adiabatic theorem the idea is to rewrite the 
integrand as a time derivative, i.e. as a commutator of H^-^^^ with an appropriately 
chosen operator A, at least up to an unavoidable error o(l). 
Let us define for n G X 

Br,{k) = i?|^(^(i/p„.(fc))Q„(fc)(i^, + fc)P„(fc) . 

Lemma 5.3. For n G T we have 

QnVfcP„ + PnVkQn = [Bn + B*. ^pcr] • 

Proof. Using the spectral decomposition and recalling 

Qn{k)VkPn{k) = ~RE„ik){Hpor{k))Qn{k){D, + k)P„{k) 

from Lemma one directly computes 
Bn{k)Hp,,{k) - Hp,,{k)Bn{k) 

= -(-ffpcr(fc) - £;„(fc))i?|„(fe)(i?pcr(fc))g„(fc)p, + fc)F„(fc) 

= -RE„(k){HpUk))Qn{k){D, + fc)P„(fc) 
= Q„(fc)VfeP„(fc). 

The lemma then follows from Pn^kQn = —{Qn^kPn)*- CH 



Lemma 5.4. 



in B{TL^TC) as e tends to zero. 



liag 



Proof. To have a concise notation in the following, expressions like W^dj'^gP„(fc) are 
understood in the sense that W^j^g acts on all fc-depending objects on its right hand 
side. We recall that W^i^g^ = PnW^Pn + QnW^Qn- Hence 

'Bn{k),Wll 

We now examine the commutators [P„(fc), W'^], [Dr^+k, W^] and [i?^ f^f,^Qn{k),W'^] 

one by one. It follows from the proof of Lemma |5.2| that [P„(fc), W^] vanishes as 
e ^ and the analogous statement for [P|, (fc)Qn(fc), W^] can be shown to hold by 
a similar argument. Thus it remains to discuss the commutator [Dx + k, W^]. For 
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i/j e H^{R'^) we compute 

W^p)iiiD^ + k) - {D, + k- p))U^){k - p) dp 

W%p)e-^p(Uijj){k -p)dp 
= e{F'm{k), 

which clearly vanishes uniformly for -0 G as e ^ 0, since F G ^fTod. vd\ 



= e 



□ 



In summary we have shown that 

>„VfcP„ + P„VfcQ„) ■F'={ Bn + Bl, H2i,^ + o(l) ) • F% 
and it remains to check 



^diag' ^ 



Lemma 5.5 

Proof. The commutator 

[Hpcr, F^] : 



in B(UH^ ,Ti,) as £ tends to zero. 



vanishes in B{H^, L?) as e — » 0. The commutator [W^J^j'"g, F'^] vanishes in B{T-L, H), 
since the commutator of P„ and Qn with F'^ are both of uniform order o(l) (in 
B{n,n)) and \W^,F^] vanishes identically. □ 

Defining 

i„ = (b„ + Bl) ■ F' , 



it follows that the integrand in (5.8) can be written as 



(Q„VfeP„ + PnVkQn) ■ F' = i^dW] + 0(1) , 

where o(l) is in the norm of B{UH^,n). (Note that for A" e B{L'^,L'^), lime^o 
A'^ = in B{L'^, L^) implies, in particular, that also linie^o A'^ = in B{H^ , L'^)). 
We are now ready for the 



Proof of Theorem S'.l. Since C^^j"g(i) : — > is bounded uniformly in t and e 
(cf. Section we obtain for the difference ( |5.8| ) of the unitary groups. 



{u^{t/e)~U^^i:^{t/e) 
(5.9) = - /* ' U^e-h - s) [i„, #"1 U'£(s) ds + o(l) . 
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Abbreviating = UH-s)U^^-i^^{s) and i„(s) = U^ii:s{-s)AnU'^Csis), we get, 

using partial integration in (5.9), 



-le 



tie 



A n" 



diag 



tie 



d ~ 



sU'it/e) I ( ) ds 



tie 



ds 



tie 



^e U^t/e) U%-s) A„ U^l^ls) ds . 



For £ ^ the first term vanishes since An is bounded and the second term vanishes, 
since tends to zero uniformly according to Lemma □ 



6. Convergence of the position operator 



In this section we will study the asymptotics of the position operator x'^{t). As 
in the case of the unitaries we have to establish that the off-diagonal contributions 
to x'^{t) vanish in the limit e — > 0. 



Proof of Theorem Let i/j e D{\x\) (1 and n el. Then 
{-'it)-Uli:,{-t/e)x^U^,i:^{t/e))i.\ 



(6.1) 
(6.2) 



< 



In order to estimate (6.1), note that we have 

rtle 



(6.3) 
and 



x%t)^p = ex^ + s / U%-s)D^U%s)^ds 



exip + £ 
exijj + e 



tie 



tie 



The last equality holds, since [Wjjj"^,^] = o(l) in B{L^), as follows immediately 



from the fact that [W^ , P„] — o{l) and [W^ , Q„] — o(l), cf. proof of Lemma |5^ 
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Hence, using ( |6.3[ ), the remaining term from ( p.l\ ) is 



(6.4) 



(6.5) 



' U-{~s/e)D, (U^is/e) ~ U'^^^^is/e)) ^ds . 



Using the fact that V and W are infinitesimally operator bounded with respect to 
-iA and that V G ff^, we get for ip{s) := U^i^^is / e)il; 

\\DlH^)\\ < H^M^) 



< 
< 



ci \\Dltl;{s)\\ + C2 U 



with ci < i and C2 < oo. Hence \\DxU^^!L{s/e)ip\\H'^ < c||'(/;||ff2 with c independent 
of s and e and we can apply Theorem |3.2| to conclude that the operator acting on 
in ( je!^ ) vanishes in B{H'^,L'^) as e ^ 0. 

We come to ( |6.5| ). Let ■(/'(s) = {U^{s/e) — U^^^{s/e))4'^ then, by Cauchy-Schwarz, 

= (V'(s),d2v-(s)) < ms)\\ \\Di^{s)\\ . 



The first factor tends to zero by Theorem 3.2 whereas the second is uniformly 
bounded by the same argument as in the treatment of (6.4) a few lines above. 
Next we rewrite (6.2) as 



with a;"^ :— QnxPn + PnxQn- This certainly vanishes as e — > if 



can be shown 



to be a bounded operator. To see this, note that in Bloch representation x acts as 
iVfc. Hence 

{UQnXPn^m = lQ,,{k)VkPn{k){U^){k) = zQ„ (fc) ( VfeF„ (fc)) (i^^) (fc) 



and thus \\QnxPn 
adjoint of Q„xP„. 



||QnVfcP„||. Finally also PnxQn is bounded, since it is the 

□ 



7. SEMICLASSICAL EQUATIONS OF MOTION FOR THE POSITION OPERATOR 

As we have shown, on the macroscopic scale the position and quasimomentum 
operators commute with the projection on isolated bands. Thus it remains to 
investigate the semiclassical limit for each isolated band separately. For this purpose 
we note that any tp £ Pn'H is of the form %l!n{k)ipn{x, k) with ipn G L'^{M*). Since 
ipn already satisfies (^^) , we have to extend the Bloch coefficients periodically. We 
determine now how H^i2^ acts on L'^{M*). We have [iJper,Pn] = and therefore 
Hpcr acts as multiplication by En{k). For 



"dC we have 



[PnW'PnU^) (k, x) 

= (27r)-'^/2 /■ w^p){ipnik),ipn{k-p))^,,^Mk-p)dpipn{k,x) 
(7.1) (M^^^>„)(fc)(^„(a;,fc). 

Thus ff^iag restricted to -P„H is unitarily equivalent to i/"^-" En{k) 
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To be able to use techniques from semiclassics we next approximate W^'" by 
the operator W^^" = W{—ieVk), where Vfc is understood with periodic boundary 
conditions on M''/r*. 



Lemma 7.1. For any n E T 

(7.2) 

in B{L^{M*)). 

Proof. By definition we have 



o{e) 



(wii^A (k) = {27t)-'/^ [ W^p)Mk - P) dp 



and therefore 
(7.3) = (27r)-'^/2 



(fc) = 



As to be shown, there exists a constant c such that 



1 ) V-'(fc ^ p)dp . 



(7.4) 



{LPn{k),ipn{k - p)) 



< c\p[ 



for Lebesgue almost all k. Therefore we conclude 



w%P)K 



?/'(fc — p)\dp 



< 



To show (7.4) note that one can chose Lpn{k) such that the map k i— > ifnik) G 
LP\M) is smooth Lebesgue almost everywhere. This is because according to Lemma 
p.l| the projections Pn{k) depend smoothly on k and hence one can locally define 
fn{k) = Pn(fc)(/5„(fco)/||Pn(fc)(/5„(fco)||. Now wc Can cover M* by finitely many open 
disjoint sets Ui such that M* \ UiUi is a set of Lebesgue measure zero and (pn{k) 
can be defined on the closure of each Ui in the way described above. One obtains a 
family ipn{k) of eigenfunctions which is smooth except at the boundaries between 
the sets, where we pick ipn{k) with an arbitrary phase. Wherever (fn{k) is smooth, 
Taylor expansion yields ipn{k — p) = (fn{k) — p ■ Vkfnik) + \p ■ M.{(pn){k')p, where 
M{ipn) denotes the Hessian and ■^p ■ M{(pn){k')p is the Lagrangian remainder. In 
view of {ipn{k),\/k^nik))L^(M) = 0, which follows from comparing ( ^.4| ) with 

(VfcP„V)(fc) = (¥'„(/=), fc))Vfc^„(fc) + (VfcV'„(fc), V'(-, k))^n{k) , 
we obtain 

'<c{k)\p\\ 



{ipn{k),ipn{k~p)) 



L2 {M) 



- 1 



Here c{k) = j \{Vn{k')Tdkidkj^n{k'))\- However, c{k) is bounded uniformly in 



k, since (fn{k) is smooth on each compact Ui. 

We define now the semiclassical Hamiltonian H^^^ 
(7.5) Hit = E,,{k) + W{-ieVk) 



□ 



acting on Lp'{M*). Then Lemma 7.1 shows that the difference iJ^'" — is of 
order 0(e) uniformly in B{L'^{M*)) and hence (cf. Section ||) the difference of the 
corresponding unitary groups approaches zero as £ — > 0. 
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then 



Corollary 7.2. Let U^^^^it) = 6"^*-^^-" and U^'''{t) = e" 

\im(U'^"{t/e)-U!'"{t/e))^0 

in B{L^{M*)). 

The semiclassical limit for U^^"'{t/e) on L'^iT*) is well studied. We refer to 
m ^ Hj. As a consequence the strong limits 

lim^ C/,^e"(-i/e) (-*eVfc) Uli^^it/e) = r„(i; fc) , 

\wUl-''\~tle)kUli''{t/e) = fc„(t;fc) 



(7.6) 
(7.7) 



e->0 



exist on iJ^(T*). r„ and fc„ act as multiplication operators and are defined as in 
(3.1) with initial conditions (r„(0), A:„(0)) = (0, A:). 

Since the restriction of ea^jjag to the n-th band subspace is unitarily equivalent to 
— ieVfc on L^(T*), we can, in view of Theorem conclude the proof of Theorem 
3.1 by showing 

Lemma 7.3. In B{L^{M*)) we have 

(7.8) Ihn {U'-'^{~t/e){-teVk)U''^Ht/s) - [/,^,"(-i/e)(-*eVfc)C/r (Ve)) = . 



Proof. The proof of (7.8) is analogous to the proof of Theorem |3.3| in Section 5, 
however, simpler. As in ( |6.3| ) we have 

U!i-i-t/e){~isVkm^^t/s) = -zeV, + s UI,-{-s) [ - zV^, i/,^] U^^is) ds 

Jo 

and 



[/^■"(-t/£)(-*£Vfc)(7^'"(V£) = 
rt/e 



C/^'"(s)ds, 



where Al^^'" := W^'" - W^/j". Now [-iVfc, = -iVkEn{k) is bounded, and 

(7.8) follows from Corollary 7.2 if we can show that [— iVfc, AVF^'"] — o(l) in 
B{L'^{M*)). Noting that (AVF'''"V)(fc) is given by (O), this can be shown by an 



argument similar to the one in Lemma 7.1 



□ 



8. Semiclassical equations of motion for general observables 
We proceed to more general semiclassical observables. First note that Theorem 



3.4 f ollow s immediately from the results obtained so far (Theorem 3.2, Corollary 7.2 
and (7.7)), since multiplication with k in Bloch representation is bounded. Hence 
we now have that 

(8.1) lim \\x^{t)i: - U~^R{t)U^\\ = 

for all -0 G RanPi r\D{\x\)nH^ and that 

(8.2) lim \\k^{t)TP - U-^K{t)U^\\ = 
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for all V' G RanPx- We next consider bounded continuous functions of x^{t) and 

Lemma 8.1. Let f £ CooO^-'^) and g G C{M*). Then for all ip e RanPj we have 

(8.3) lim II {fix^t))-Vl-'fiRit))U)^P\\=0 

and 

(8.4) lim II {g{k^t)) - U-'g{K{t))U) ^|| = . 

Proof. We will sketch the proof for x^it). First note that R{t) := U-'^R{t)U is 
a bounded self-adjoint operator and commutes with Pi. Hence the sets D± := 
{R{t) ± i)(RanPi n £>(|a::|) n H"^) are dense in Pj (Since R and x'^ are vectors of 
operators in M'', note that this and the following statements hold component wise). 
For ij} £ D± we have 

(8.5) [(x^{t) ± - {R{t) ± i)-^] ip = {x'{t) ± tr\Rit) - x%t))ip 



for If = {R{t)±i)-^ij G RanFinL'dxDni?^. Thus, by Theorem |3.2|, ( p.5|) strongly 
approaches zero as e — * and, since D± are dense in Pi, (a;^(t) ± i)^^ strongly 
approach {R{t) ± i)^^ on Px. 

Using the fact that polynomials in {xj ± i)"^ , j — 1, . . . , d, are dense in C^aiW^) 
one concludes that the convergence [t) R{t) on RanPj in the "strong resolvent 
sense" implies 

lnn||(/(a;-(t))-/(P(t)))V||=0 

for aU / e Coo(K'^) and -0 G RanPr (cf. Theorem VIII.20 in ||||). However, by the 
functional calculus for self-adjoint operators we have f{U^^R{t)U) = U^^ f{R{t))U 



and (8^) f ollow s 



Clearly (8.4) follows analogously. □ 



Proof of Theorem 8.5. Let a G 0(0). Referring again to the general Stone- Weier- 
straB theorem we can uniformly approximate a(a;,^) by a sum of products, i.e. 
ai^.O = EZo'^Mx)g,iO with /, G C^iR''), g, S C(M*), E|a,| < oo and 

we are 



supjgpj^gjfd^gAf* \fii^)9iiO\ < Hence in order to prove Theorem 3.5 
left to show that for arbitrary / G Coo(]R'^) and g G C{M*) we have 

(8.6) (t) - U-'f{R{t))g{K{t))U 

strongly on RanPj. To see this recall the so called product rule for quantum 
observables (cf. |2^). It states, in particular, that for two symbols A,B<E 0(0) 

lini II {{AB)^-^ - A^-^B^^") = . 
Applied to our case this yields 

ifi^)9i0)'^^' (t) - ifixr-'aiO'^^') (t) - fix^t))gik^t)) . 
Finally, since / and g are bounded. Lemma |8.l| implies (^.6[) and thus Theorem 



3.5. □ 
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9. Band crossings 

We proved the semiclassical limit for isolated bands only. In principle, there are 
two distinct mechanisms of how this assumption could be violated. First of all a 
band could be isolated but have a constant multiplicity larger than one. This occurs, 
e.g., for the Dirac equation where because of spin the electron and positron bands 
are both two-fold degenerate. A systematic study is only recent ^ and leads to 
a matrix valued symplectic structure for the semiclassical dynamics. For periodic 
potentials degeneracies are the exception. They form a real analytic subvariety of 
the Bloch variety B = {{k, X) G R'^ x R \ 3f G L^{M) : H^.{k)_f = A/} and have 



a dimension at least one less than the dimension of B \17, Thus points of 

band crossings have a fc-Lebesgue measure zero. From the study of band structures 
in solids one knows that band crossings indeed occur. Thus it is of interest to 
understand the extra complications coming from band crossings. 

There are two types of band crossings. The first one is removable through a 
proper analytic continuation of the bands. In a way, removable band crossings 
correspond to a wrong choice of the fundamental domain. E.g. for V — we 
may artificially introduce a lattice F. The bands touch then at the boundary of 
M*. Upon analytic continuation we recover the single band Ei{k) = with 
M* =W^. In one dimension all band crossings can be removed |^. Thus, with the 
adjustment discussed, our result fully covers the case d = 1. For d > 2 generically 
band crossings cannot be removed. 

It is then of great physical interest to understand how a wave packet tunnels 
into a neighboring band through points of degeneracy (or almost degeneracy) . For a 
careful asymptotic analysis in particular model systems we refer to the monumental 
work of G. Hagedorn Gerard considers a model system with two bands in 

two dimensions, i.e., the role of —iA + V^ is taken by ^ ^fc ^ ' Reinvestigates 

the semiclassical limit and proves that the particle may tunnel to the other band 
with a probability which depends on how well the initial wave packet is concentrated 
near a semiclassical orbit hitting the singularity. 
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